IFT/9/1994 



EQUATION OF STATE FOR PARTICLES ARISING AT THE UNIVERSE 
AT GRAND UNIFICATION ENERGIES 



I. Dymnikova 

N. Copernicus Astronomical Center 
'^f ' Bartycka 18, PL-00-716 Warsaw, Poland 

ON ' 

ON ; N. Copernicus Foundation for Polish Astronomy 

Al. Ujazdowskie 4, PL-00-478 Warsaw, Poland 

>V 

1 M. Krawczyk 

Institute of Theoretical Physics, Warsaw University 
>^ Hoza 69, PL-00-681 Warsaw, Poland 

m 
m 
in 
o 

Abstract 

ON 

First postinflationary stage of the Universe evolution is considered in more detail. It 
is shown that heavy particles with mass Mh ~ M-gvt arising at the Universe at phase 
transitions at Grand Unification Energies behave like ideal quantum degenerate Bose gas. 
The equation of state for both scalar and gauge bosons is presented including the coupling 
constant and vacuum expectation value at Equt- One possible way is proposed to connect 
cosmological observational data with parameters of Grand Unified Theories. 
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The equation of state is necessary instrument for constructing the cosmological mod- 
els. Putting it into the Friedman equations, one gets cosmological solution describing a 
certain stage of the evolution of the Universe (see, e.g. 1 ). In the case of an inhomogeneous 
situation, the equations are more complicated, but procedure is on essence the same (for 
review see 2 ). 

On the other hand, the equation of state can also serve as useful tool for testing physics 
of very high energies. Calculating cosmological model with an equation of state at Grand 
Unification energy Equt ~ 10 15 GeV one can hope to get an opportunity to extract some 
information concerning GUT physics from cosmological observational data 3 . 
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In the cosmological context the equation of state has the form 



P = «£ 7 > (!) 

where p is pressure and e energy density of matter. The case a = corresponds to the 
dust equation of state, the case 7 = 1 and a = 1/3 describes the ultrarelativistic particles, 
and the case 7 = 1 and a = — 1 describes vacuum 4 (or false vacuum 5 ). 
The vacuum, or false vacuum, or inflationary equation of state 

p=-e (2) 

has been proposed about thirty years ago. In 1965, considering various possible equations of 
state at superhigh densities, Sakharov 6 has proposed (2) as equation of state for superdense 
matter. At the same time Gliner 4 has proposed (2) as the equation of state for vacuum 
with nonzero energy density, because corresponding stress-energy tensor 

T a p = zgotp (3) 

has an infinite set of comoving reference frames and no distinguished one. Gliner has 
suggested that at superhigh densities a superdense matter undergoes a transition into 
a state of superdense vacuum. In 1970 he also suggested, that such a state could be 
initial state for the expanding Universe with negative pressure as the intrinsic reason for 
expansion 7 . Corresponding cosmological model has been calculated in 1975 8 : Because 
initial vacuum nonsingular state (2) is homogeneous and isotropic, and expansion is starting 
from the causally connected region confined by the de Sitter event horizon ao, here denoted 
by r : 

^r 4 

the homogeneity and isotropy of the Universe have been guaranteed in obvious way. 

In 1981 Guth has found that eq.(2) can arise in Grand Unified Theories as the equa- 
tion of state for scalar field in metastable state called false vacuum 5 , and now it is broadly 
used in inflationary scenarios (for review see 9 ' 10 ). Since 1981 many various and sophisti- 
cated physical mechanisms have been proposed driving the inflation and providing further 
Universe reheating (for review 9 ' 10 ). 

In the present paper we consider the first postinflationary stage. The state (2) is 
unstable in both classical 7 and quantum 1 ' context. As a result of phase transitions 9 , 
heavy nonrelativistic particles arise with mass Mh ~ Mqut 9 ' 10 - Note, that light particles 
can not live in the causally connected region (4), which at the energy scale ~ 10 15 GeV is 
of a size 10 -25 cm. 

The first phase transition can be related to the spontaneous global symmetry breaking 
(because arrow of time is expected to appear just then), hence at first scalar bosons emerge 
at the Universe. So, we are starting with and focusing upon scalar bosons. On the other 
hand, all our estimates below are on essence valid for both scalar and gauge bosons, and 
we present equation of state for both of them. 
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First heavy particles emerging at the Universe can not be at rest, at least because 
they are quantum particles, and they can undergo collisions. As we do not expect that 
they appear with large momenta, the CM energy of their collisions \/s is close to the sum 
of rest energies of colliding particles y/s ~ 2Eh (the threshold of the reaction) . Therefore 
the corresponding interactions are mainly due to elastic 2^2 body scatterings. The 
production of one extra particle will be possible at much higher CM energy, greater then 

To describe our heavy particles by the equation of state, we should consider them as 
gas. To learn what kind of gas they form, we should determine characteristic radius of 
interaction for them (depending on the coupling constant) and compare it to geometrical 
radius, which for quantum particles is defined by the Compton wavelength: 



M H c 



Xcompt ~ n— ~ 10 29 cm. (5) 



Corresponding effective geometrical cross-section is defined by 

= K\ 2 Cornpt . (6) 

Assuming that first heavy particles are scalars with the selfinteraction, we can estimate 
their characteristic radius of interaction by the mass of particles mediating the interaction. 
Since these particles are the same scalars with masses M H , the radiation of such a particle 
for 2 — > 2 body scatterings is possible only due to the nonconservation of the energy 

AE ~ M H c 2 . (7) 

Then, at the threshold, the cross section for the interaction mediated by heavy particles 
and depending on the coupling constant gn is given by 

(Jint ~ olh 1 /M% ~ a H 2 &o, (8) 

where 

a H = g 2 H /4n. (9) 
Now let us define the characteristic radius of interaction by the relation 

Vint = nr 2 int . (10) 

Note , that Uj n t <C a® , if an is small, thus the characteristic radius of interaction in 
considered case is even smaller than the Compton wavelength. 

Heavy particles emerge in the vacuum with the energy density 

M GUT , 4 77 , 3 , . 

Pvac = PGUT = PPlanckiYJ ) ~ 10 9/071. (11) 

M Planck 

The radius of causally connected region defined by the de Sitter formula (4) is equal to 

r ~ 10" 25 cm. (12) 
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Hence we can estimate the average momentum of the particles as 



\p\ ~ h/r ~ 10 12 GeV (13) 

and average kinetic energy as 

E kin = \p\ 2 /2M H ~ h 2 /2M H r 2 . 

For the ratio we have 



E kin /E H ~ % 2 /{r 2 Q M 2 H c 2 ) ~ (Acom P t/ro) 2 ~ 10" 8 (14). 

So, we can consider our heavy particles as dust, at the first approximation. To tell more 
about them we should make sure that we have ensemble of them to use statistics. This is 
just moment when we need some cosmological model of transition, and now we know that 
we should organize transition from vacuum (2) into the dust with pressure p = (up to 
10~ 8 according to (14)). 



We describe the phase transition using phenomenological cosmological model. Follow- 
ing 8 ' 15 , let us assume that during a transition the equation of state has the form 

(15) 

At e = eo we have vacuum equation of state (2). At e = e± we have dust equation of 
state p = 0. A parameter < f3 < 1 characterizes the rate of transition, for example, if a 
transition is driven by scalar field , (3 characterizes the rate of rolling down a scalar field 
potentials with small or nonexisting barrier between false and true vacuum. For simplicity 
let us choose (3 = 1/2 and the case of the spatially-flat model (k = 0). 

Integrating the Friedmann equations with the transitional equation of state (15) we 
get for the scale factor a the expression 

ct 

a = cio ex.p{A sin — — }, (16) 
Aa 

where 

l-o ^ 1 '^ ( 17 ) 

and for the energy density the expression 

£ = £ (l-sin 2 ^-). (18) 

Ado 

Using condition of conservation of the rest mass at the end of the transition, afei = 
a today £ today, 1 , we get e\ ~ (l/65)ff0) ^i ~ 10 _6 cm and duration of the transition t\ ~ 
2 x 10~ 34 sec, adopting p to day ~ 10~ 29 g/cm s . 
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Because we organized a transition into the dust, the whole rest mass (Mi ~ 2 x 10 57 s) 
is contained in our heavy particles with Mh ~ Mqut ~ 2 x 10~ 9 g. So, number of particles 
is Ni ~ 10 66 . The volume at the end of transition is V± ~ 4 x 10 _18 cm 3 , hence density of 
particles is 

n = ^ „ 10 84 cm- 3 . (19) 

Using the Friedmann equations and formulae (16)-(18), we can check the total energy 
balance during the transition and find that only about 1/4 of vacuum energy involved in 
transition has been spent. The substantial amount remains available for further expansion 
and reheating. Hence, the Universe is still mainly vacuum dominated, and particles have 
at their disposal to interact the causally connected region (12). How many particles are in 
this region? 

Npart.in caus.conn.reg . ^caus.conn.reg * ^ 

Vcaus.conn.reg. = K 2 rl (de Sitter volume). 

So, 

Npart.in caus.conn.reg. ~ 10 • (20) 

Hence, we have an ensemble and opportunity to use statistics. 



According to (19), the mean distance between particles is 

r av ~ 10" 28 cm. (21) 

The characteristic radius of interaction defined by (10), r^ n t ~ XcomptctH is much smaller 
than r av . Since the characteristic cross section (8) is smaller than geometric cross-section 
(6), the free path is determined by 

I = — 3 x l(T 27 cm. (22) 

As a result our particles satisfy two general criteria for ideal gas: 

nr int ^1 ~~ * First criterion for ideal gas 
Tint <^ I —> Second criterion for ideal gas 

Note, that the free path I is much less than the radius of causally connected region: 

/ < r . (23) 

Hence, in caussaly connected region particles can collide and produce some temperature 
and pressure. 

Next question is what kind of ideal gas we have. The average momentum is given by 
(13), then the de Broglie wavelength can be estimated as 

XdeBroglie = V^T ~ 10" 25 Cm (23) 
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Because XdeBroglie ^ r av, our ideal gas is essentially quantum one. Indeed, the above 
condition means that number of particles N is much more than the number of quantum 
states at their disposal 



( XdeBroglie , 3 _ N h_ _ N 

1 r av 1 ~V p* - {Vp 3 )/n s » 



With | p | 2 ~ MnkTn where k = 1.38 x 10 16 er5r/if - the Boltzman constant, the condition 
(24) takes more frequently used form 16 

n h 3 

V{M H kT H f/2 >>L W 

So, we have ideal quantum Bose gas. 

Its temperature of degeneration is determined by 16 

Ti 2 

kT deg = 3.31-27= n 2 '\ (26) 

Here g st = 2s + 1 is statistical weight, with s = for scalar and s = 1 for vector bosons. 

For Mh ~ 2 x 10 — 9 (7 and n ~ 10 84 cm -3 we obtain /cTd e£/ ~ 10 11 er(7. For our gas with 
average momentum defined by (13), kTn ~ E^in ~ 3 x 10 erg. Hence, its temperature 
Th Tdegi and we can surely describe our gas by the equation of state for the degenerate 
Bose gas 16 with the pressure given by 

M 3/2 

P = 0.08blg st ^(kT H ) 5 / 2 (27) 

According to the standard approach (with selfinteraction A</> 4 , A = g 2 H ) to generating 
masses of heavy scalars via the spontaneous global symmetry breaking 13 ' 14 ? we have 

M H ~ v l9H , (28) 

where vj is a vacuum expectation value. To reveal characteristic dependence of pressure 
on the coupling and the vacuum expctation value, we can put i) average kinetic energy for 
kTii and ii) express the quantity ro in terms of Mpi and Mqut- Then the equation (27) 
takes the form 

P = const C ^-^-(^fM" GUT . (29) 
n vig H Mpi 



We expect that Vj ~ Mqut, so 



U gn Mpi 



We see that in the equation of state for quantum ideal degenerate Bose gas the coupling 
constant is appearing. Because, as it can be easily estimated, pressure is much less than 
energy density, the process of arising of particles looks like evaporation of Bose-condensate. 
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It seems to be natural to expect that next in the earliest evolution of the Univers 
some gauge interaction appears with gauge bosons X, Y and their gauge coupling constant 
g. There are estimation of the value of this gauge coupling at the energy scale 10 15 GeV, 
based on the accelerator's data at low energies 17 . It gives for the gauge coupling constant 

14,17 

a = g 2 /ATT ~ 1/40- 1/25. (31) 

The spontaneous local symmetry breaking is expected to give a mass of these gauge 
bosons of the order 13,14 

M x ~ gv H . (31) 

Since both effects are at similar scale ~ Equt, and this is just scale of our analysis, 
we can put vi ~ vu and gn ~ g- (Note, that in SUSY case vj = vu and gu = g-) Hence 
the equation of state for the gauge bosons will be the same as for scalar ones (see (27) or 
(29,30)) with appropriate change in constant. 

The results presented here do not depend in a crucial way whether the scale Mqut 
is equal to 10 1 5 or (what is preferable by SUSY models) 10 1 6 GeV. There is also possible 
to have both phase transitions I and II at the same time. 

We can calculate cosmological model including equation of state (29) depending on 
GUT parameters 3 . That way we could get opportunity to extract some information 
concerning GUT physics from observational cosmological data, and compare it with accel- 
erator's data. It would provide some way toward astronomical testing GUT models. 
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